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Abstract 

In this paper, we present the results of Monte Carlo simulations of 7-ray and cosmic ray 
proton induced extensive air showers as detected by the TACTIC atmospheric Cherenkov 
imaging telescope for optimizing its trigger field of view and topological trigger generation 
scheme. The simulation study has been carried out at several zenith angles. The topological 
trigger generation uses a coincidence of 2 or 3 nearest neighbour pixels for producing an 
event trigger. The results of this study suggest that a trigger field of 1 1 x 1 1 pixels (~ 3.4° 
X 3.4°) is quite optimum for achieving maximum effective collection area for 7-rays from 
a point source. With regard to optimization of topological trigger generation, it is found 
that both 2 or 3 nearest neighbour pixels yield nearly similar results up to a zenith angle 
of 25° with a threshold energy of ~ 1.5 TeV for 7-rays. Beyond zenith angle of 25°, the 
results suggest that a 2-pixel nearest neighbour trigger should be preferred. Comparison of 
the simulated integral rates has also been made with corresponding measured values for 
validating the predictions of the Monte Carlo simulations, especially the effective collec- 
tion area, so that energy spectra of sources (or flux upper limits in case of no detection) 
can be determined reliably. Reasonably good matching of the measured trigger rates (on 
the basis of ~ 207 hours of data collected with the telescope in NN-2 and NN-3 trigger 
configurations) with that obtained from simulations reassures that the procedure followed 
by us in estimating the threshold energy and detection rates is quite reliable. 
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1 Introduction 



Detection of cosmic 7-rays in the very high energy range ( ~ 100 GeV-50 TeV) us- 
ing ground-based atmospheric Cherenkov technique [1-2] allows a detailed study 
of the nonthermal universe and probes many phenomena of fundamental astrophys- 
ical interest including the acceleration and emission mechanisms in galactic and 
extragalactic sources, physics of pulsars and supernova remnants, study of extra- 
galactic diffuse background etc. The recent success of the field has mainly resulted 
from the development of the Cherenkov imaging telescopes which allows efficient 
separation of photon induced showers from the hadron background. The power of 
this technique has been convincingly demonstrated by several modem atmospheric 
Cherenkov telescopes ( e.g. MAGIC [3], VERITAS [4], HESS [5] and CANGA- 
ROO [6]) and there is now a general consensus that stereoscopic arrays provide the 
most promising avenue for future observations. 

A typical imaging atmospheric Cherenkov telescope consists of a light collector 
and a close-packed array of fast photomultiplier tubes ( also called the imaging 
camera with individual tubes as its pixels). The imaging camera records the spa- 
tial distribution of photons in the image (called the Cherenkov image) and by ex- 
ploiting the differences in the images produced by 7-ray and cosmic-ray initiated 
showers, caused by the physical processes responsible for their development in the 
atmosphere, it becomes possible to effectively segregate the two event species on a 
statistical basis with a high degree of accuracy. Simulation studies of atmospheric 
electromagnetic (7-ray induced) and hadronic ( proton induced) cascades [7,8] have 
shown that there exist subtle differences in the image features ( shape, size and ori- 
entation) of these event types. The image of an off-axis 7-ray shower is relatively 
more compact and is elliptical in shape, with the major axis of the ellipse oriented 
towards the camera centre (i.e. the source position). In the case of off-axis cosmic 
ray shower, on the contrary, the image is comparatively bigger and irregular, be- 
sides having no preferred orientation (because of isotropic nature of cosmic rays). 
Remarkably indeed, these predictions have been thoroughly matched by observa- 
tions of various groups with great success and it is primarily because of the success 
of these experiments that the field of ground-based 7-ray astronomy today boasts a 
catalog of more than a hundred 7-ray sources [9]. 

In an atmospheric Cherenkov imaging telescope, there are several crucial param- 
eters which need to be fine-tuned for maximizing the performance of the system. 
Assuming that we are observing a point 7-ray source with the telescope axis di- 
rectly pointing towards it, these are : (a) the point spread function of the light col- 
lector, (b) the pixel size of the imaging camera, (c) trigger field of view and (d) 
trigger generation scheme. The optimum trigger field of view of an atmospheric 
Cherenkov imaging telescope needs to be chosen in such a manner that the ratio 
of 7-rays to the cosmic ray background is maximum. Thus, for an imaging camera 
with a large field of view, it is generally advantageous to exclude the outer pixels 
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from the trigger generation. These pixels only help to properly reconstruct the im- 
age parameters and do not increase the 7-ray rate. The trigger in an atmospheric 
Cherenkov imaging telescope is generated by demanding a coincidence between a 
number of PMTs, with the condition that pixels participating in the trigger should 
have a signal exceeding some threshold qo. The conventionally used multi-fold ( 2 
to 4) coincidence approach necessitates the operation of the telescope at a relatively 
higher threshold energy in order to keep the chance coincidence rate at a reason- 
able value. An improved trigger generation scheme involving two or three nearest 
neighbour pixels can be used instead for lowering the threshold energy of the tele- 
scope. The main aim of this work is to use Monte Carlo simulations for optimizing 
the above mentioned parameters for the TACTIC telescope. 



2 TACTIC imaging telescope 

The TACTIC (TeV Atmospheric Cherenkov Telescope with Imaging Camera) 7- 
ray telescope [10] has been in operation at Mt. Abu ( 24.6° N, 72.7° E, 1300m asl), 
India, for the last several years to study TeV gamma ray emission from celestial 
sources. The telescope uses a tessellated light-collector of area ~ 9.5m^ which is 
capable of tracking a celestial source across the sky. The telescope deploys a 349- 
pixel imaging camera, with a uniform pixel size of ~ 0.31" and a ~ 5.9° x 5.9° 
field-of-view, to record atmospheric Cherenkov events produced by an incoming 
cosmic-ray particle or a 7-ray photon with an energy above ~lTeV. The TACTIC 
light-collector uses 34 front-face aluminium-coated, glass spherical mirrors of 60 
cm diameter each with a focal length ~ 400cm and reflection coefficient >80% at 
a wavelength of ~ 400nm. The point-spread function has a FWHM of ~ 0.185° 
(= 12.5mm) and D90 ~ 0.34" (=22. 8mm). Here, D90 is defined as the diameter of 
circle, concentric with the centroid of the image, within which 90% of reflected rays 
lie. The value of D90 ~ 0.29° (=19. 3mm), predicted on the basis of the simulation 
for an incidence angle of 0°, matches reasonably well with the measured value 
mentioned above. Other details regarding the ray-tracing simulation procedure and 
comparison of the measured point-spread function of the TACTIC light collector 
with the simulated performance of ideal Davies-Cotton and paraboloid designs are 
discussed in [1 1]. The photographs of the TACTIC imaging telescope and its back- 
end signal processing electronics are shown in Fig.l 

The imaging camera uses 19mm diameter photomultiplier tubes (ETL-9083 UVB). 
The bialkali photocathode has a maximum quantum efficiency of ~ 27% and for 
this PMT the maximum quantum efficiency is usually reached in the wavelength 
range ~ 340 -360 nm. In order to ensure better light collection efficiency and re- 
duction in the background light falling on the photomultipliers we use Compound 
Parabolic Concentrator (CPC) as light guides on the photomultipliers. The CPCs 
used in the TACTIC imaging camera have the following specifications : entry aper- 
ture ~21.0 mm; exit aperture ~15.0 mm; and height ~17.6 mm. The light col- 
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Fig. 1. (a) Photograph of the 349-pixel TACTIC imaging telescope (b) Photograph of back- 
-end signal processing electronics used in the telescope. 

lection efficiency of the CPCs is estimated to be of ~65% by using ray tracing 
simulations and includes both the geometrical collection efficiency and the reflec- 
tivity of the surface. The light collection efficiency mentioned above refers to those 
photons which undergo reflection after falling on the inside surface of the CPCs. 
The data acquisition and control system of the telescope [12] has been designed 
around a network of PCs running the QNX (version 4.25) real-time operating sys- 
tem. The triggered events are digitized by CAM AC based 12-bit Charge to Digital 
Converters (CDC) which have a full scale range of 600 pC. The relative gain of the 
photomultiplier tubes is monitored by repeatedly flashing a blue LED, placed at a 
distance of ~1.5m from the camera. Whenever the single channel rate of any two 
or more pixels in the trigger region goes outside the preset operational band, it is 
automatically restored to within the prescribed range by appropriately adjusting the 
high voltage of the pixels [13]. 

The imaging camera uses a programmable topological trigger [14] which can pick 
up events with a variety of trigger configurations. As the trigger scheme is not 
hard wired, a number of coincidence trigger options ( e.g Nearest Neighbour Pairs, 
Nearest Neighbour Non-coUinear Triplets and Nearest Neighbour Non-CoUinear 
Quadruplets) can be generated under software control. The number of pixels par- 
ticipating in the trigger can be varied from 7x7 pixels to 15x16 pixels. Apart from 
generating the prompt trigger with a coincidence gate width of ~18ns, the trigger 
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generator has a provision for producing a chance coincidence output based on C2 
combinations from various groups of closely spaced 12 channels. It is worth men- 
tioning here that there are several reasons which compel us to use a coincidence 
gate width of ~18 ns instead of an optimal gate width < 5 ns. The reasons are the 
following: (a) variation in the transit time of the PMTs because of being operated 
at different high voltage values, (b) variation in the propagation delay of 50m long 
RG58 cables because of procurement of the cables in different batches, (c) degrada- 
tion in the rise time of the PMT pulses when the signal is transmitted through a 50m 
long RG58 cable and the subsequent 'time- walk' of the leading edge discriminators 
which is dependent on the amplitude of the PMT output signal and (d) variation in 
the access time of static RAMS ( Toshiba make TC55B417) used in the trigger 
generator and difference of ~ 4ns in the cascading and the non-cascading routes 
[14]. Feasibility of using programmable digital delay generators for fine tuning the 
signal propagation time in all the channels is also being investigated so that the 
coincidence gate can be reduced to < 10ns. Other details regarding the design, im- 
plementation and performance evaluation of the programmable topological trigger 
generator for the 349-pixel imaging camera of the TACTIC telescope are discussed 
in [14]. 

A user friendly in-house developed tracking system software provides an indepen- 
dent movement for the zenithal and azimuthal axes so that a matching between the 
telescope pointing direction and the source direction is obtained with an accuracy 
better than ±2 arc-min. The tracking accuracy of the telescope is also checked on 
a regular basis with so called 'point runs', where a reasonably bright star, having a 
declination close to that of the candidate 7-ray source is tracked continuously for 
about 5 hours. The point run calibration data (corrected zenith and azimuth angle 
of the telescope when the star image is centered) are then incorporated in the tele- 
scope drive system software so that appropriate corrections can be applied directly 
in real time while tracking a candidate 7-ray source. Other relevant details of the 
instrumentation aspects of the telescope are discussed [10]. 



3 Importance of using neighbouring pixels for trigger generation and TAC- 
TIC trigger generator design 

The conventional approach of generating an event trigger in atmospheric Cherenkov 
imaging telescopes is to use a majority type coincidence with trigger multiplicity of 
2 to 4. The optimum value of the trigger multiplicity is decided by several factors 
like observatory altitude, primary energy and zenith angle dependent intensity dis- 
tribution of the Cherenkov photons in the image plane and the pixel size. However, 
in order to keep the chance coincidence rate at a reasonable value (usually few per- 
cent of prompt coincidence rate), the use of a majority type coincidence approach 
( i.e. any m present out of total n detector channels, denoted by any-m,n here on- 
wards ), in generating the trigger, leads to a relatively higher threshold energy. The 
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reason for this being the presence of "C^ factor in m-fold chance coincidence rate 
(Jany-m,n) formula which is given by 

-1 -1 

f any— m.n ~ C(xny—m,n'^Rpixel\QO) 'Tr ^ 77 \T -^pixeliQo) c 

m!(n — my. 

where Cany-m,n= "'Cm= n\/{m\{n — m)\) is the number of all possible m-pixel 
combinations, Rpixeiilo) is the shot-noise generated single channel rate passing the 
discriminator set at a threshold of qo (in units of photoelectrons), m is the number 
of detector channels which fire coincidently within a coincidence resolving time 
Tc, n is the total number of detector channels which can participate in producing 
the trigger. It has been assumed in the above formula that the contribution from 
higher multiplicity triggers (i.e. due to coincidence of > m channels) is negligible 
as compared to the m-fold chance coincidence rate and this can be satisfied easily 
as long as Rpixei{Qo)Tc <<!• Keeping in view the compact nature of 7-ray images, 
an elegant way to reduce the chance coincidence rate is to limit the number of 
m-fold combinations by using topological constraints in the trigger (such as strict 
adjacency of triggering pixels or demanding that the pixels should lie within camera 
sectors). If CNN-rn,n IS thc number of m-fold nearest neighbour combinations the 
corresponding chance coincidence rate (fNN-m,n) is given by : 

fNN-m,n ~ C]siN-m,n1^Rpixel{(lo)^'^^ ^ (2) 

Since, the number of nearest neighbour m-fold combinations ( CNN-m,n) is always 
less than the corresponding number of combinations in a simple majority type coin- 
cidence ( Cany-m,n), thcrc is a significant reduction in the chance coincidence rate 
when nearest neighbour m-fold trigger generation is used. Thus, use of a nearest 
neighbour m-fold based trigger generation allows the individual pixels to operate at 
a relatively lower threshold which eventually leads to lowering the threshold energy 
of the telescope. 

An improved trigger generation scheme, based on the above concept was proposed 
by us [15] for the Whipple system and a reduction in threshold energy by about 
20% is achieved by the group when a 2 adjacent /331 trigger condition instead of 
the conventional ^^^6*2 trigger condition is used [16]. This has resulted in ~40% 
increase in the number of 7-rays detected from the Crab Nebula [16]. The useful- 
ness of the nearest-neighbour trigger generation scheme in achieving the lowest 
possible energy threshold has been demonstrated by other telescopes also includ- 
ing VERITAS, HESS and MAGIC [17,18]. The TACTIC imaging telescope also 
employs a programmable topological trigger generator with 2( or 3) nearest neigh- 
bour pixels participating in the trigger and the main aim of this work is to optimize 
the overall trigger generation scheme, both with respect to trigger multiplicity (i.e. 
NN-2 versus NN-3) as well as number of pixels participating in the trigger. Thus, 
the work involves comparing the performance of NN-2 and NN-3 trigger schemes 
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with number of pixels participating in the trigger varying from 49 to 240. 

The number of 2-fold and 3-fold combinations for majority type and Nearest- 
Neighbour coincidence for different numbers of pixels participating in the trigger 
is given in Table. 1 These combinations have been worked out by considering the 

Table 1 



The number of 2-fold and 3-fold combinations for majority and Nearest-Neighbour coin- 
cidence trigger schemes. 



No. of Trigger 
pixels (n) 


Field of 
View 


Number of Possible combinations 


2-fold 


3-fold 






Cany—2,n 


CNN-2,n 


Cany—3,n 


CNN-3,n 


7x7 (=49) 


2.17° x 2.17° 


1176 


156 


18424 


624 


9x9 (=81) 


2.79° x 2.79° 


3240 


272 


85320 


1152 


11x11 (= 121) 


3.41° X 3.41° 


7260 


420 


287980 


1840 


13x13 (= 169) 


4.03° X 4.03° 


14196 


600 


790244 


2688 


16x15 (=240) 


4.96° X 4.65° 


28680 


869 


2275280 


3968 



actual layout of the pixels used in the TACTIC imaging camera ( Fig.2a). Different 
trigger regions with 49,81,121,169 and 240 pixels considered in this work are also 
indicated in this figure. Defining r/^ „ = CNN-m,n/ Cany-m,n as the reduction in the 
chance coincidence rate. Fig. 2b shows the variation of ?72,n and rj^^n as a function 
of n. With r]2,n ^ 4.822 n'^-^^^ and 773,^ ~ 44.182 n-^-^"^^', as seen'from Fig.2b, it 
is obvious that using a nearest neighbour trigger approach offers an elegant way to 
reduce the chance coincidence rate. Allowing only certain m-pixel combinations to 
participate in the trigger also helps in removing a pixels (or certain pixels) with ex- 
cessive single pixel rate from the trigger without making hardware changes. A very 
common reason for the excessive single pixel rate is due to presence of relatively 
bright stars in the field of view of the telescope ( e.g. a star like ^-Tauri when Crab 
nebula is being tracked by a telescope). 

Regarding the design aspects of the TACTIC trigger generator, while a detailed de- 
scription of the topological trigger used in the TACTIC telescope can be found in 
[14], we will only present here some of its key design features. The programmable 
topological trigger used in the TACTIC telescope is based on 16k x 4 bit fast static 
RAMs and has a coincidence resolving time of ~18ns. As the trigger scheme is 
not hard wired, a number of coincidence trigger options ( e.g NN-2, NN-3 etc.) can 
be generated under software control. Given the rectangular layout of the TACTIC 
imaging camera pixels we use a maximum of 240 inner pixels for trigger gen- 
eration. Since each 16k x 4 bit RAM can accommodate only 16 pixels, we have 
divided the trigger pixels into 20 cells of 12 pixels each in the form of a 3x4 ma- 
trix. The outputs from each of the four independently operating trigger generator 
modules are then collated in an Event Handler module which generates the 22ns 
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Fig. 2. (a) Layout of the 349-pixels of the TACTIC imaging camera. Different trigger 
regions with 49, 81, 121, 169 and 240 pixels considered in this work are also indicated in 
the figure, (b) Variation of r]2,n and 773 „ as a function of number of pixels in the trigger 
region. 

duration gate pulse and interrupts the Data Acquisition System (DAS) for reading 
the charge ADC data from all the 349 pixels. 



4 Behaviour of single pixel rate and effect of afterpulsing in PMTs 



It is a very well known fact that the trigger threshold of an atmospheric Cherenkov 
telescope is dictated by afterpulsing in PMTs instead of fluctuations in the number 
of photoelectrons due to the light of night sky (LONS) background alone [19,20]. 
The main aim here is to determine the single pixel threshold in a realistic manner 
so that the same can be used in Monte Carlo simulations for reliably predicting the 
telescope performance ( e.g. trigger efficiency as a function of energy, differential 
rates etc.) when it responds to 7-ray and cosmic ray showers. The phenomenon of 
afterpulsing in PMTs leads to generation of large amplitude pulses (well in excess 
of 10 photoelectrons) at a non-negligible rate. Thus, it is the single pixel rate and 
the resulting maximum permissible chance coincidence rate which decides the min- 
imal discriminator threshold and ultimately the energy threshold of the telescope. 
While the phenomenon of afterpulsing (also known as ion feedback) has been in- 
vestigated thoroughly by a number of workers, we will only discuss here, for the 
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sake of completeness, some important consequences of this effect. Afterpulsing 
is a result of ionization of a molecule or an atom of the residual gas inside the 
PMT (or adsorption at dynode surface material) by photoelectron and subsequent 
release of many photoelectrons from the photocathode when the positively charged 
ion hits the photocathode. The relevance of afterpulses depends on how the PMTs 
are used and they can largely be inconsequential if the main pulse can be gated. 
Since the PMTs used in atmospheric Cherenkov imaging cameras are operated in 
a self-trigger mode (i.e. the PMTs are live all the time), afterpulses pose a serious 
problem. 

Before presenting the measured behaviour of single pixel rate as a function of 
threshold level, it is worth estimating the single pixel rate ( Rpixei) as a function 
of the threshold level on the basis of Poissonian fluctuations in the number of pho- 
toelectrons. Assuming that the individual LONS photons (and hence the photoelec- 
trons) are uncorrelated and the discriminator threshold is set at qo photoelectrons, 
RpixeiiQo) can be estimated by the following formula : 

Rpixeiiqo) ~ Rpe ex.-p{-Rpers) (3) 



where is the characteristic shaping time of the single photoelectron response of 
the PMT ( approximated as a rectangular pulse); Rpe is the single photoelectron 
rate at which photoelectrons are generated; Rpixei is trigger rate such that, within 
a time r^, there are qq photoelectrons. It has been also assumed in equation 3 that 
RpixeiiQo) ^ Rpe and contribution from > (go) electrons is negligible to the trigger 
rate. The single photoelectron rate (Rpe), for a given light intensity, can be obtained 
by monitoring the average PMT current (/„) and the PMT gain (Gpmt) and then 
using the following formula: 

Rpe = (4) 



where e is the electron charge. An estimate of the threshold level ({go)), on the basis 
of the equation 3, can now be readily worked out for a particular telescope system. 
Using a light of night sky background intensity ~ 2 x 10^^ photons m'^s^sr~^ in 
the wavelength range ~320-650nm for the Mount Abu observatory and other rele- 
vant hardware specifications of the TACTIC telescope ( e.g., light collection area, 
pixel solid angle, quantum efficiency of the PMT etc. ) we find that average number 
of photoelectrons per 10ns time interval is ~0.20. With a single photoelectron rate 
of ~20 MHz, on using equation 3, it turns out that the PMT trigger rate should be 
~ 1 kHz and ~2.3 x 10^^ Hz at threshold levels of 5 and 10 photoelectrons, respec- 
tively. Actual measurements of the PMT rate as a function of threshold, however, 
yield much higher values because of afterpulsing. Realizing that the discriminator 
threshold is a very crucial input required for performing Monte Carlo simulations 
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in a realistic manner it becomes necessary to take actual trigger rate measurements. 
Accordingly, we have also carried out these measurements for several PMTs and 
the results of this study are discussed below. 

For determining the afterpulsing behaviour, laboratory measurements were carried 
out for several PMTs by simulating the LONS with a LED. The intensity of the 
LED was set in such a manner, that the resulting PMT anode current is equivalent to 
the value which the PMT would produce during actual observations. A summary of 
the results obtained from these measurements is shown in Fig. 3. The integral count 
rate as a function of threshold for 1 PMT using a LONS background ~ 2.0 x 10^^ 
photons m~'^s~^sr~^ is shown in Fig.3a ( marked as 10.3 //A). In order to simulate 




Fig. 3. (a) Variation of trigger rate as a function of discriminator threshold for 3 representa- 
tive values of Ught of night sky background values, (b) Differential count as a function of the 
discriminator threshold for determining the single photoelectron spectrum, (c) Variation of 
trigger rate as a function of discriminator threshold for 3 different gain calibrated PMTs. (d) 
Measured amplitude distribution of shot-noise fluctuations in the output of a PMT exposed 
to a steady Ught source (Uke the light of night sky). 

the behaviour of integral rate when the telescope is observing a relatively brighter 

or fainter region in the sky, we also show 2 more representative curves ( marked as 
15.4 ^A and 6.7 ^A) in this figure . The corresponding LONS background values 
for these two curves has been chosen to be ~ 3.0 x 10^^ photons m^'^s^^sr'^ and 
~ 1.3 X 10^^ photons m~'^s~^sr~^, respectively. It is important to mention here 
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that, during actual observations, we operate the PMTs at a lower gain and the typ- 
ical values of anode currents are < 3 jiA for most of PMTs. Operation of PMT at 
higher gain is only followed during laboratory measurements so that the PMT can 
be calibrated easily. The calibration of the discriminator threshold in terms of pho- 
toelectrons (top X-axis scale of Fig.3a ) has been done by determining the gain of 
the PMT using standard single photoelectron detection method. The measurement 
for this purpose was taken under very low illumination level and the differential 
count rate as a function of discriminator threshold is shown in Fig.Bb. The integral 
count rate measurements for 3 additional gain calibrated PMTs is shown in Fig. 3c 
and it is evident from this figure that the afterpulsing rate can vary significantly for 
the same type of PMT. Although a variation of only ~± 25% was reported by us 
in our earlier work [14] on the basis of about 5 PMTs from the same batch, de- 
tailed measurements performed, over last several years indicate the variation in the 
afterpulsing rate, especially for PMTs which are from different batches, is signifi- 
cantly larger than this value. Factors like variation in the vacuum quality, purity of 
the dynode material used and the differences in the high voltage needed to reach 
the same gain can be some of the reasons responsible for this variation. The differ- 
ences in the afterpulsing rate of PMTs also depends on the threshold level where 
they are compared. For a large sample of PMTs, procured in different batches, the 
differences in the afterpulsing rate can be even by a factor of 10 in the threshold 
level range of ---'10-20 pes. Taking guidance from [21], a function of the type Rap ~ 
RpeXa exp{—qob) ( where Rap is the afterpulsing rate ; a and b are fit parameters) 
has also been fit to the afterpulsing rate of the 3 PMTs ( Fig. 3c) for (go) > 5 pe. The 
best fit values of b obtained are ~ 0.29, 0.49 and 0.56. A typical value of b~0.5 has 
also been reported for the same PMT type by the HEGRA group [21]. 

Referring back to Fig. 3c, one can readily use this plot to determine the discrimina- 
tor threshold for a predetermined value of PMT trigger rate. For example, taking 
the case of NN-2 trigger with CMN-m,n = 420 (i.e. 121 pixels in the trigger region), 
Tc ~18ns and an additional condition that fNN-m,n < 0-04 Hz ( ~ 2% of expected 
cosmic ray rate of ~2Hz ); by using equation 2 we find that Rpixei{%) should be in 
the range of 50 Hz to 100 Hz with corresponding values of go ~22pe. On repeating 
the above exercise for NN-3 trigger with CNN~m,n = 1840 it is seen that RpixeiiQo) 
should be around 3000 Hz with a corresponding threshold go ~ 13pe . Keeping in 
view that several key telescope parameters (e.g. effective collection area, threshold 
energy, differential and integral rates) depend on the value go, the above information 
about go provides a very important input for carrying out Monte Carlo simulation 
studies in a realistic manner. The values of go, for the TACTIC telescope with NN-2 
and NN-3 trigger configurations have also been marked in Fig. 3c . Although one 
can determine the value of go by performing Monte Carlo simulations at several 
trial go values and eventually choose the one for which the integral rate estimates 
match the experimental results, following this approach requires intensive compu- 
tational resources. While we have also followed this approach but the main aim of 
the exercise was only for the purposes of fine tuning the value of go and best results 
were obtained for the following trigger threshold values : go ~22pe for NN-2 trig- 
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ger and go ~14pe for NN-3 trigger. We will thus only use these values of qo for rest 
of the work described in the paper. 

In order to carry out realistic simulations for the TACTIC telescope, the simu- 
lated data-bases of the photoelectron content of each pixel in a Cherenkov image 
also needs to include the shot-noise contribution from night sky background light. 
Fig. 3d shows the measured differential pulse amplitude distribution of the shot- 
noise fluctuations in the output of a PMT. The gate pulse to the charge to digital 
converter unit for this measurement was generated by using a pulse generator. Us- 
ing an appropriate probability density function for this noise profile curve, we have 
then superimposed the effect of shot-noise on the Cherenkov light generated pho- 
toelectron content of each pixel in the camera. 



5 Generation of Monte Carlo simulated data base for TACTIC telescope 
using CORSIKA code 

For the present studies, simulations were carried out using the CORSIKA air- 
shower code (Version 5.6211) [22] with Cherenkov option. The CORSIKA sim- 
ulation code simulates interactions of nuclei, hadrons, muons, electrons, photons 
as well as decays of unstable secondaries in the atmosphere. It utilises the EGS4 
code for the electromagnetic component of the air shower simulation and dual par- 
ton model for the simulation of hadronic interactions. The simulations have been 
carried out for Mt.Abu observatory altitude of ~1300 m with appropriate values of 
35.86 f/T and 26.6 f/T respectively for the horizontal and the vertical components 
of the terrestial magnetic field. 



In order to keep the data-base size and computation time within manageable limits, 
we have performed multiple sampling of a single shower at several core distances. 
A total of 104 detector locations (matrix of 52 x 2 elements), with each detector hav- 
ing a size of 4m x 4m, has been used so that a TACTlC-like telescope can sample 
Cherenkov light at the core distances of upto 245 m along one cartesian axis. The 
position of the shower core has been fixed at the geometrical center of the 52x2 
array. The centre-to-centre inter-element spacing is 10m along the x-direction and 
6 m along the y-direction. We have also taken care of the zenith angle dependence 
on the size of each detector element while performing the simulation at various 
zenith angles. Again, in order to keep the data-base size within manageable lim- 
its, the Cherenkov photon wavelength band chosen for generating the database is 
300-450 nm, and the Cherenkov photon bunch size used is varied from 1 to 30 (de- 
pending on the energy and arrival direction of the primary particle). The simulated 
data-base for 7-ray showers uses about ~ 38035 showers in the energy range 0.2- 
20 TeV with an impact parameter up to 250m. These showers have been generated 
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Table 2 



The number of gamma and proton progenitor showers used in simulation for different 
values of primary energies at various zenith angles. 



Particle Type 
Gamma / Proton 


Number of showers at various zenith angles 












Energy (TeV) 


5° 


15° 


25° 


35° 


45° 


0.2/0.4 


1000 / 1000 

1 \J\J\J 1 1 \J\J\J 


1000 / 1000 

1 \J\J\J 1 1 \J\J\J 


0/0 


0/0 


0/0 


0.3 / 0.6 


1 000 / 900 


1 000 / 900 


1 000 / 1 000 

1 \J\J\J 1 1 \J\J\J 


0/0 


0/0 


0.4 / 0.8 


1 000 / 800 


1 000 / 800 


1 000 /1 000 

1 \J\J\J 1 1 \J\J\J 


/o 


0/0 


0.55/1.1 


1 000 / 750 

J. \J\J\J 1 1 0\J 


1 000 / 700 

J. \J\J\J 1 1 XJXJ 


1 000 / 1 000 


1 000 / 1 000 

A-\J\J\J 1 1-\J\J\J 


0/0 


0.7/1.4 


900 / 700 

y\j\j 1 1 \J\J 


900 / 600 


900 /900 


900 / 900 


900 / 1 000 


1.0/2.0 


800 / 650 


800 / 550 


800 / 800 


800 / 800 


800 / 900 


1.3/2.6 


660 / 600 


700 / 500 


700 / 700 

/ \J\J 1 1 \J\J 


700 / 700 

/ \J\J 1 1 \J\J 


700 / 800 

/ \J\J 1 \J\J\J 


1.8/3.6 


400 / 550 


400 / 450 


500 / 600 


60f) / 600 


6f)0 / 700 


2.5 / 5.0 


325 / 500 


350 / 400 


400 / 500 


500 / 500 


500 / 600 


3.3/6.6 


300 / 450 


325 / 350 


350 / 450 


400 / 450 


400 / 550 


4.5 / 9.0 


275/400 


300 / 300 


325/425 


350 / 425 


350 / 500 


6.0/12.0 


250 / 350 


275 / 250 


300 / 400 


325 / 400 


325 / 450 


8.0/16.0 


225 / 300 


250 / 200 


275 / 375 


300/375 


300/400 


11.0/22.0 


200/275 


225 / 175 


250/350 


275 /350 


275 / 350 


15.0/30.0 


175 / 250 


200/150 


225 / 325 


250 /325 


250/300 


20.0/40.0 


150/225 


175 / 125 


200 / 300 


225 / 300 


225 / 250 


Total 


8660 / 8700 


8900 / 7450 


8225/9125 


6625/7125 


5625 / 6800 



at 5 different zenith angles ( i.e. 9= 5°, 15°, 25°, 35° and 45°). Since it is known 
that the contribution of cosmic ray protons is highest to the total event rate of at- 
mospheric Cherenkov telescopes [23] with other nuclei, including He, contributing 
up to ~30% of the total event rate, we have used only proton primaries in our sim- 
ulation study. Taking guidance from the work [24], where a similar approach has 
been followed in the simulations, we have scaled up the event rate calculated for 
protons by a factor of ~ 1.5 in order to get a total cosmic ray rate. Accordingly, 
a data-base for proton showers uses about ~ 39200 showers in the energy range 
0.4-40 TeV and the incidence angle of these showers is simulated by randomizing 
the arrival direction of the primary in a field of view of 6.6° x 6.6 ° around the 
pointing direction of the telescope. The number of 7-ray and proton showers simu- 
lated for different primary energies and at various zenith angles are listed in Table 
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2. It is important to mention here that the number of gamma-ray showers as well as 
the number of proton showers have not been generated according to a power law 
distribution. The main reason for this is to avoid the problem of limited statistics 
at higher energies on account of too few showers. However, appropriate 7-ray and 
proton spectra have been folded in the results while determining differential and 
integral rates. Similarly, the total number of events(or images) considered for the 
analysis work is sampled in an off-line manner as per ~ r dr scaling, where r is the 
distance of the telescope from the center of the array. The total number of images 
considered for the present study is ~ 1.9 million for 7-rays and ~ 2.0 million for 
protons. 

The CORSIKA code generates a Cherenkov file containing seven parameters for 
characterizing the Cherenkov photon (or photon bunch) received with a detector 
element area (4mx4m in our case). These parameters are bunch number, x and 
y co-ordinates, two direction cosines, height of production and arrival time at the 
ground. This CORSIKA code which we have used does not consider wavelength 
dependent atmospheric absorption of the Cherenkov photon (bunch). Accordingly, 
atmospheric absorption of the Cherenkov photon-bunch as well as its subsequant 
history, is taken care of by a supplementary BACUP-code, developed in house. 
Essential details related of this code are described below. 

The parameters of each photon-bunch, produced in the wavelength range, A ~ 300- 
450 nm by simulation code CORSIKA and managing to the reach 52x2 detector 
array, are written in the Cherenkov file. The BACUP-code uses this information to 
obtain corresponding parameters of photons expected in the overall wavelength re- 
gion A ~ 200-600 nm. Keeping in view the variations in the atmospheric extinction 
coefficient at shorter and longer wavelength, the step-size used is 5 nm, 10 nm, 20 
nm and 50nm for the ranges A = 200-270 nm, 270-280 nm, 280-400 nm and A > 
400 nm respectively. The number of Cherenkov photons expected in various wave- 
length intervals is extrapolated using ~ A~^ dependence of the Cherenkov light 
emission spectrum. Each Cherenkov photon is traced back to the point of emission 
and the A- dependent atmospheric extinction is considered for the wavelength range 
A = 200-270 nm, and A > 270nm, using measured data and using data derived from 
[25] respectively. Each photon-bunch, after correcting for atmospheric absorption, 
is reflected from the mirror-facet and ray-traced to the corresponding position in 
the focal plane of the telescope. Appropriate care is taken to include the reflectivity 
of mirrors, transmission of the light guides (CPC) and quantum efficiency of the 
phototubes in the present simulation studies. A modest value of ~82% for the pho- 
toelectron collection efficiency has been used in teh simulation. The obstruction 
encountered due to the telescope mechanical structure by the incident and reflected 
photon-bunches, during their prorogation, is also duly considered. The data-base 
of resulting photo-electrons (pe) in the imaging camera pixels, are subjected to 
noise injection and image cleaning for further analysis. Representative examples of 
Cherenkov images recorded at a zentih of angle 25° is shown in Fig. 4. The grey 
shade color scheme followed, for representing number of photoelectrons recorded 
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X (deg) X (deg) 

\ ^^^^^ ^^^ 

0.00 0.60 1.21 1.81 2.41 3.02 3.62 4.22 4.83 5.43 

ln( photoelectrons) 

Fig. 4. (a) Cherenkov image of a 7-ray shower with energy ~4.5 TeV and a core distance 
of ~ 85m. (b) Cherenkov image of a proton shower with energy ~9.0 TeV, core distance 
of ^ 85m and arrival angle ~2.87°. 

in a pixel, using equal In(photoelectrons) intervals, is also shown at the bottom of 
this figure. The Cherenkov image of a 7-ray shower, shown in Fig.4a, has a primary 
energy of ~4.5 TeV and a core distance of ~ 85m. The Cherenkov image, shown in 
Fig.4b, for an off-axis proton primary has the following parameters : energy ~9.0 
TeV ,arrival direction ~2.87° and core distance of ~ 85m. 



6 Formalism for calculating effective collection rate, threshold energy and 
detection rates 

The trigger rate of an atmospheric Cherenkov imaging telescope and the resulting 
statistics of the showers is primarily decided by two basic parameters; effective 
collection area and the energy threshold. For a given configuration of the telescope, 
these parameters depend on the pixel trigger threshold, the trigger multiplicity and 
the trigger field of view. The energy threshold of the telescope primarily depends on 
the geometrical area of the light collector and the photon to electron conversion ef- 
ficiency. The later being determined by the mirror reflectivity, collection efficiency 
of the light guides, the spectral response of the PMTs and other possible losses of 
the Cherenkov light on the way to the PMTs. 

The effective collection area of an atmospheric Cherenkov imaging telescope is 
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governed by the lateral and the angular distribution of the Cherenkov light gener- 
ated by the primary particle (i.e. 7-ray or a cosmic-ray). For a 7-ray, coming from 
a point source, the effective collection area ( A^{E)) is given by 

00 

A^{E) = 27ry" P^{E,r)rdr (5) 




where P-yiE, r) is the detection probability (or trigger efficiency) for a 7-ray with 
primary energy E and impact parameter r. By the trigger efficiency -Ry(-E', r), we 
imply the fraction of the simulated Cherenkov events which can be detected by the 
telescope placed at a given distance r from 7- ray or proton- generated EAS core. 
Since, for the isotropic cosmic-ray background the detection probability Pcr{E, r, O) 
depends also on the angle between the direction of the primary particle and the tele- 
scope axis (6), the effective collection area ( Acr{E)) is given by 

2n 00 

Acr{E) = '271 J J PcR{E,r,e)rdrdfl (6) 




where dil = 2n sin 6 and the integration over 6 is included in dVt. It is important to 
mention that the effective collection area for cosmic rays, given by equation 6, has 
units of rri^ sr. Using the effective collection area derived from above equations, the 

differential rate can be found by multiplying the effective area with the correspond- 
ing differential spectra of 7-rays and cosmic -rays ; i.e dR^J dE = A^ (E) x dF^/dE 
and dRcR/dE = Acr{E) x dFcR/ {dEdVt). For 7-rays we have used the Crab 
Nebula differential spectrum [26] 

= 6.0 X lQ-^E-''-^^-'''''''''<^)m-^s-^TeV-^ (7) 



For the cosmic-ray background, we have used the following spectrum [24] 

( ^^^A = Qm&E-^'m-^-hr-^TeV-^ (8) 
\dEdVt 



As already mentioned in section 5, we only use proton showers for cosmic -ray 
simulations. The reason for this being the fact that for a given energy, the show- 
ers induced by heavier nuclei develop at substantially larger heights in the atmo- 
sphere with a less intense Cherenkov light pool. In addition, the fraction of en- 
ergy channeled into electromagnetic subshowers, responsible for the production of 
Cherenkov photons decreases with increasing nucleon number. A combination of 
all these result in a larger but less intense Cherenkov light pool which increases 
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the threshold energy of the telescope for heavier cosmic ray particles. The integral 
trigger rate for 7-rays (Rj) and cosmic-rays (Rcr) is given by 



00 



dK 




dE 



A^{E) dE 



(9) 







00 



dFcR 
dEdn 



Rcr — 



I 



Acr{E) dE 



(10) 







It is important to mention here that we have discretized the integrals ( equations 5 
and 6) in our simulations with a step size of 10m in core distance and a step size of 
~0.1° for the arrival direction. The largest value of core distance used in the sim- 
ulation is ~ 245m ( imposed by using an array of 490 m length in the CORSIKA 
simulations) and this restriction yields reasonably accurate values for primary en- 
ergy range considered in the present work. Significant deviations in calculation of 
the effective area are expected at higher primary energies ( > 40 TeV), but the rela- 
tive contribution of such particles is only marginal. 



7 Optimization of the trigger field of view 

While pixel size and the total field of view are the basic parameters in the design of 
a camera for an imaging telescope, it is also important to optimize the trigger field 
of view so that, 7-ray showers with an impact parameter of up to 150m, are able to 
trigger the telescope. Keeping in view that the detection rate of showers initiated by 
isotropically distributed cosmic rays is proportional to the solid angle, using a trig- 
ger field of view beyond the optimum value only reduces the ratio of 7-ray trigger 
rate to cosmic -ray trigger rate. The optimum trigger field of view is decided by the 
position of the image centroid in the camera and is determined by the height of the 
shower maximum above the observatory level and the impact parameter. Knowing 
that, for 7-ray showers, the angular shift of the image from the camera center can 
be estimated quite reliably using a simple toy model of the electromagnetic cascade 
development, it is worth having a look first at these estimates before proceeding to 
detail Monte Carlo simulations. In the toy model, one assumes that most of the 
Cherenkov light photons are coming directly from the shower maximum. If H^ax 
is the height of shower maximum (with Xmax as the corresponding atmospheric 
depth) where the number of secondary electrons in the shower is maximum, the 
expression for Xmax, using classical cascade theory is given by : 



(11) 
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where is the primary energy of the 7-ray; Ec is the critical energy (Ec ~ 76 



MeV) and Xo is the radiation length in the atmosphere (xo ~ 36.6 gcm~ ). The 
angular displacement of the image centroid from the camera center depends on the 
impact parameter and the line of sight distance to the shower maximum and is given 
by the following expression [2]: 



tan 



bcos9 



bsinO + /loln[1030/(x^aa;COs6')] - Hobs 



(12) 



where ho is the scale height of the atmosphere (ho ^ 8.5 km ) ; Ho^s is the obser- 
vatory altitude; b is the impact parameter and 9 is the zenith angle. The variation of 
angular shift of the image centroid as a function of impact parameter for different 
7-ray energies is shown in Fig.5. On examining this figure, it is clearly seen that 




Fig. 5. Angular displacement of the image centroid from the camera center as a function of 
core distance for various zenith angles ( 9= 0°, 20° and 40°) and primary 7-ray. (a) Ej=l 
TeV (b) E^=3 TeV (c) £^^,=5 TeV and £;^=10 TeV. 

an imaging camera with a trigger field of view of ~ 1.5° (half-angle) is able to de- 
tect 7-ray showers of energies < 5TeV up to impact parameter values of ~150m. 
It is also seen in the figure that the angular shift of the image centroid from the 
camera center decreases when zenith angle is increased. Furthermore, if the main 
concern is to maximize the ratio of 7-ray trigger rate to cosmic ray trigger rate, 
further increase in the trigger field of view does not lead to an increase in the 7 -ray 
trigger rate because of the fact that the detection of 7-ray showers beyond ~ 150m 
is limited due to the reduced density of the Cherenkov photons. 
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For optimizing the trigger field of view on the basis of Monte Carlo simulations we 
have used simulated data generated at two zenith angles (viz., 15" and 35") with 
trigger regions having 49, 81, 121, 169 and 240 pixels. The corresponding fields 
of view of these trigger regions are ~ 2.17° x 2.17", - 2.79" x 2.79°, ~ 3.41° x 
3.41", ~ 4.03" X 4.03" and ~ 4.65" x 4.96°. Both NN-2 as well as NN-3 trigger 
generation schemes, with a single pixel threshold ( qo) of ~ 14pe and ~ 22pe, 
respectively, have been considered in the simulation study. The trigger efficiency 
as a function of 7-ray energy, for showers falling with a core distance of ~250m 
is shown in Fig.6. Referring first to Fig. 6a and Fig. 6b ( i.e efficiency curves at a 
zenith angle of 15" for NN-2 and NN-3 trigger modes), it is evident that at > 4.0 
TeV, the trigger efficiency increases when the number of pixels in the trigger region 
is increased from 49 to 121 . It is also seen that the trigger efficiency remains almost 
constant when the number of trigger pixels is increased further from 121 to 240. A 
similar behaviour is also evident in Fig.6c and Fig.6d when a zenith angle of 35° is 
used for studying the trigger efficiency. Therefore, these results suggest that using 




Fig. 6. Trigger efficiency at a function of energy of the primary 7-rays with trigger regions 
having 49, 81, 121, 169 and 240 pixels, (a) NN-2 trigger with 6'=15", (b) NN-3 trigger with 
9=15" (c) NN-2 trigger with 9=35° and (d) NN-3 trigger with 9=35°. 

121 pixels ( 11x11 matrix), having a field of view of ~ 3.41" x 3.41", is quite 
optimum for producing triggers in response to 7-rays with an impact parameter up 
to ~ 150m. 

We have also estimated the rate loss for cosmic rays when the largest value of 
core distance used in the simulations is restricted to ~245m. This restriction arises 
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because of using an array of 490m length in the CORSIKA simulations. In this 
study we have included showers with an impact parameter up to 300m and the 
behaviour of the cosmic ray rate as a function of core distance is shown in Fig.7(a- 
d). Referring first to Fig. 7a and Fig. 7b ( i.e cosmic ray rate at a zenith angle of 
15" for NN-2 and NN-3 trigger modes), one can easily see that, for a trigger field 
of view with 240 pixels, there is only a marginal increase in the cosmic ray rate ( 
~2.5%) when the core distance increases from ~245m to ~300m. It is also seen 
in the figure that corresponding increase in the cosmic ray rate is negligible for a 
trigger field of view of with 49 pixels. A similar behaviour is also seen in Fig.7c 
and Fig.7d when a zenith angle of 35" is used. 



(a) NN-2| 6 = ^5°^—^ 
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Fig. 7. Cosmic ray rate as a function of core distance with trigger region having 49, 121 
and 240 pixels, (a) NN-2 trigger with 9=15°, (b) NN-3 trigger with 9=15° (c) NN-2 trigger 
with 9=35° and (d) NN-3 trigger with ^=35". 



The performance of the imaging camera with 49, 8 1 , 1 2 1 , 1 69 and 240 pixels in the 
trigger region has also been studied by calculating the integral 7-ray and cosmic -ray 
trigger rate and the results of this study are shown in Fig. 8. The 7-ray trigger rates 
shown in Fig. 8a and Fig. 8b have been calculated by using the Crab nebula energy 
spectrum. The corresponding trigger rates for cosmic-ray protons as a function of 
trigger field of view are shown in Fig.Sc and Fig.Sd. Strong dependence of the 7- 
ray trigger rate as a function of trigger field of view is only seen when 49 pixels are 
used for trigger generation. A weak dependence of 7-ray trigger rate, for all other 
configurations with > 121 pixels in the trigger region, is also evident from these 
figures. With trigger rates of ~ 40 events h~^ for NN-2 trigger and ~ 30 events 
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Fig. 8. (a & b) Gamma-ray trigger rate as a function of trigger field of view for NN-2 and 
NN-3 trigger schemes, (c & d) Cosmic-ray trigger rate as a function of trigger field of view 
for NN-2 and NN-3 trigger schemes. 



h"^ for NN-3 trigger at a zenith angle of 15", it is also seen from these figures that 
NN-2 trigger scheme provides better detection efficiency at higher zenith angles 
as compared to the NN-3 trigger scheme. Linear dependence of the proton trigger 
rate on the trigger field of view seen in Fig. 8c and Fig.Sd, in accordance with the 
isotropic behaviour of cosmic rays, also reassures that the trigger field optimization 
results are reliable. Therefore it can be safely concluded that a trigger region using 
innermost 121 pixels in the TACTIC imaging camera provides an optimal field of 
view for detecting 7-rays from a point source. 

It is important to mention here that the actual field of view of the camera should 
be more than the trigger field so that truncation of the Cherenkov images can be 
avoided [21]. The shapes of Cherenkov images, produced by 7-rays, get otherwise 
distorted and the effect is especially important for the LENGTH parameter (de- 
fined as the rms angular size and is related to the longitudinal development of the 
shower). The limited field of view also introduces large errors in the determination 
of other image parameters and in estimating the energy of the primary 7-rays. It 
is for this reason that imaging camera of the TACTIC telescope uses 225 to 349 
pixels. Further work related to realizing the potential of the full camera with 349 
pixels, especially using Wobble mode of observation is still in progress. 
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8 Optimization of topological trigger generation scheme on the basis of esti 
mated threshold energy and detection rates 



After having validated that a trigger field of view with 121 pixel is quite opti- 
mum for detecting the 7-rays from a point source, we now proceed to performance 
evaluation of the TACTIC telescope on the basis of detection rates of 7-rays and 
cosmic-rays. Optimization of the topological trigger condition by considering 2 
and 3 fold nearest neighbour trigger schemes (i.e. NN-2 versus NN-3) is another 
aspect which is discussed in this section. Collection areas for 7-ray and cosmic-ray 
induced showers can be derived from Monte Carlo simulations by evaluating the 
energy dependent trigger efficiency {P^{E, r) and Pcr{E, r, 0)). Using the equa- 
tions 5 and 6, along with a single pixel threshold of ~ 22pe for the NN-2 trigger 
and ~ 14pe for the NN-3, the Monte Carlo evaluated effective collection areas, as a 
function of the progenitor particle-energy, are shown in Fig.9 for 3-different zenith 
angles. The results have been obtained by using the inner-most 121 pixels (11x11 
matrix) for trigger generation. To facilitate an easier comparison of effective ar- 




Fig. 9. Effective collection areas for 7-ray and proton induced showers, at various inci- 
dent angles, with NN-2 and NN-3 trigger schemes, (a) 7-rays with NN-2 trigger scheme 
(b)7-rays with NN-3 trigger scheme, (c) Protons with NN-2 trigger scheme.(d) Protons 
with NN-3 trigger scheme. 

eas, we have divided the effective collection area for protons by Cltrig, where Vltrig 
denotes the trigger field of view solid angle (~ 3.54 x 10"^sr). The energy depen- 
dence of the effective collection area on the primary energy is clearly evident in 
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this figure and reflects the shape of the lateral distribution of Chrenkov photons on 
the ground. Compared with the relatively steep lateral distribution of the cosmic- 
ray showers, the lateral distribution of Cherenkov light produced by 7-ray showers 
has a rather flat lateral distribution up to ~120m core distance ( a more precise 
value is dependent on the observatory altitude). Using derived effective collection 
areas shown in Fig. 9, the corresponding differential rates, for 7-rays and protons is 
shown in Fig.lO(a-d). The differential rate curves, for the two event species, have 




Energy (TeV) Energy (TeV) 

Fig. 10. Estimated differential rates for 7-ray and proton induced showers, at various inci- 
dent angles, with NN-2 and NN-3 trigger schemes ; (a) 7-rays with NN-2 trigger scheme 
(b)7-rays with NN-3 trigger scheme; (c) Protons with NN-2 trigger scheme ;(d) Protons 
with NN-3 trigger scheme. 

been obtained by using Crab Nebula spectrum for the 7-rays (equation 7 ) and pro- 
ton spectrum ( equation 8). Following the standard definition of effective threshold 
of the telescope i.e., value of the primary energy at which the differential rate peaks, 
the resulting energy thresholds of the TACTIC telescope are presented in Table 3. 
Results regarding integral rates, calculated by using equations 9 and 10, are shown 
in Fig. 11 and also given in Table 3. It is important to mention here that we have 
scaled up the calculated event rate for protons by a factor of ~ 1.5 in order to esti- 
mate the total cosmic ray rate ( as discussed in section 5). The main purpose of this 
scaling is to account for contributions from other nuclei, including He. 

On comparing the performance of NN-2 and NN-3 trigger configurations, one can 
notice the following : (i) In the case of NN-3 trigger, there is a marginal increase 
in the 7-ray threshold energy when zenith angle is increased from 5° to 25°. (ii) 
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Fig. 11. Estimated integral rates for 7-ray and cosmic ray induced showers, at various in- 
cident angles, with NN-2 and NN-3 trigger schemes, (a) 7-rays with NN-2 trigger scheme 
(b)7-rays with NN-3 trigger scheme, (c) Cosmic-rays with NN-2 trigger scheme.(d) Cos- 
mic-rays with NN-3 trigger scheme. 

Table 3 

Threshold energy for 7-rays (Crab-nebula spectrum) and protons for the TACTIC tele- 
scope, with 121 pixels in the trigger region and using NN-2 and NN-3 trigger configura- 
tions. The integral detection rates for 7-rays and cosmic-rays (CR) are also given in the 
table. A scaUng factor of ~ 1 .5 ( as discussed in section 5) has also been used to estimate 
the total cosmic-ray rate. 



Zenith angle 
(in degree) 


Energy Threshold (TeV) 


Integral Rate 


NN-2 


NN-3 


NN-2 


NN-3 


7 


Proton 


7 


Proton 


7(h-i) 


CR(s-i) 


7 (h-') 


CR (s -1) 


5 


1.3 


2.6 


1.3 


2.6 


40 


2.0 


40 


2.4 


15 


1.3 


3.6 


1.8 


3.6 


36 


1.9 


34 


2.1 


25 


1.3 


3.6 


1.8 


3.6 


50 


1.6 


34 


1.8 


35 


1.8 


3.6 


2.5 


3.6 


36 


1.3 


22 


1.4 


45 


3.3 


9.0 


6.0 


9.0 


18 


1.0 


11 


1.0 



Beyond 6>25" the 7-ray threshold energy increases more rapidly for NN-3 trigger 
(~6.0TeV at 45°) as compared to the NN-2 trigger. In the case of NN-2 trigger the 
7-ray threshold energy changes from l.STeV at ^=35° to ~3.3TeV at ^=45°. The 
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resulting integral trigger rate for 7-rays as a function of zenith angle is also seen to 
drop more rapidly in NN-3 trigger mode compared to the NN-2 trigger configura- 
tion. The reason for this being the fact that at larger zenith angles the Cherenkov 
images produced by 7-rays tend to be more compact as a result of increase in the 
line of sight distance to the shower maximum. Using a NN-2 trigger methodol- 
ogy is thus a better option as against the NN-3 trigger configuration when a source 
needs to be observed at 9> 30°. Advantage of using a NN-2 trigger configura- 
tion, for the HEGRA telescope with a pixel size of ~ 0.25°, has also been pointed 
out in [21] while comparing the performance of various nearest neighbour trigger 
schemes with different multiplicities. 



9 Comparison with experimental data 

The 349-pixel TACTIC imaging telescope, has been in operation at Mt.Abu, India 
since 2001 and has so far detected 7-ray emission from the Crab Nebula, Mrk 
421 and Mrk 501. First successful detection of TeV 7-rays from the Crab Nebula 
was seen in the data collected by the telescope for ~ 41.5h between January 19 - 
February 23, 2001 . The data was taken by using a NN-3 topological trigger with the 
innermost 240 pixels (15x16 matrix) of the camera participating in generating the 
event-trigger. A statistically significant excess of ~6.3(t was seen in the data with 
~447ih71 7-ray like events. During the last few years, regular observations were 
taken on a number of potential 7-ray sources (viz., 1ES2344-I-514, PSR 0355-1-54, 
ON 231, H1426, Mrk-421, Mrk-501 etc.) and the results of these observations have 
been presented in [27-30]. A major portion of this data has been collected with 
121 pixels (11 X 11 matrix) participating in the trigger with either NN-2 or NN-3 
trigger configuration and using a total of 225 pixels. At present, the telescope has 
a 5a sensitivity of detecting Crab Nebula in 25 hours of observation time. This 
sensitivity figure has been verified several times over last 5 years and consolidated 
results on the Crab nebula, for ~310.5 h of observation time taken over a period of 
3 observing seasons, is in preparation.. 

For validating the Monte Carlo simulation results for NN-2 trigger mode we have 
used data on the Crab Nebula for ~101.44 h between November 10, 2005 - January 
30, 2006. In this data sample, the number of 7-ray events obtained after applying 
the Dynamic Supercuts are found to be ~(928±100) with a statistical significance 
of ~9.40(j [31]. Similarly, for validating the simulation results for NN-3 trigger 
mode we have used data on the Crab Nebula for ~ 105. 16 h between November 
03, 2007 - March 8, 2008. After applying the Dynamic Supercuts to this data sam- 
ple, the number of 7-ray events are found to be ~(968±90) with a statistical sig- 
nificance of ~ 11.1(7. Data collected in the NN-2 trigger mode, during November 
10, 2005 - January 30, 2006, has also been used by us previously for checking 
the validity of the ANN-based energy estimation procedure and for determining 
the energy spectrum of the Crab Nebula in the energy range 1-16 TeV [31]. The 
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7-ray differential spectrum obtained after using appropriate values of effective col- 
lection area and 7-ray acceptance efficiency (along with their energy and zenith 
angle dependence) was found to be follow a power law {d^/dE = foE~^) with 
/o = (3.12 ± 0.48) X IQ-'^^cm-^s-'^TeV-^ and T = 2.69 ± 0.14 [31]. While the 
details of the standard Dynamic Supercuts procedure used by us to separate 7-ray 
like images from the background cosmic rays can be found in [27,31], in this work, 
we will only concentrate on the cosmic ray data collected by the telescope for val- 
idating the Monte Carlo simulation estimates regarding effective collection area, 
threshold energy and integral rates. 

Starting with the analysis of the NN-2 trigger mode data collected for ~ 101. 44 h; 
we have first divided the total data into 4 spells, where each spell corresponds to 
one lunation period. Without correcting for inter-pixel gain variation and without 
applying other data selection cuts. Fig. 12a shows measured values of cosmic-ray 
rate as a function of zenith angle. The zenith angle plotted in this figure has been 
multiplied with the sign of hour angle so that pre-upper transit and post-upper tran- 
sit rates can be distinguished easily. Each point in this plot represents a 10 minute 
observation run and the event rate has been calculated by dividing the number of 
events recorded in 600s. The behaviour of the measured values of cosmic-ray rate 
as a function of zenith angle for the NN-3 trigger mode is shown in Fig. 12b. In 
this case, the total data has been divided into 5 spells. The corresponding Monte 
Cario estimates (solid line), with RcR{d) ~ 2.012{cosdy-^^^H z for NN-2 trigger 
and Rcr{0) ~ 2.330(cos^^)^'^^°iJ2; for NN-3 trigger, are also shown in Fig. 12a and 
Fig. 12b, respectively. Although the Monte Carlo estimates are found to describe the 
zenith angle dependence of the measured rates reasonably well with an accuracy of 
± ~20% for most of the 10 minute runs, presence of a systematic uncertainty in the 
experimental data can not be ruled out on account of sky condition variations, both 
on daily basis as well as from one lunation to other. Differences in parameters like 
single pixel trigger threshold, overall photoelectron conversion factor etc. are also 
some of the additional factors because of which there can be a possible mismatch 
between simulated trigger rates and the corresponding experimental values. While 
work on understanding the telescope systematics is still in progress, we have used 
a value ~ 15% for systematic error in the measured cosmic-ray rate for further 
analysis. The justification for this is based on the guidance from [32] where the au- 
thor has reported a variation of ~ 15% — 20% in the Chrenkov photon density for 
mid-latitude sites on account of seasonal changes. Although, one does not expect a 
significant variation in atmospheric extinction on a nightly basis for a good astro- 
nomical site, significant variation in atmospheric extinction on a nightly basis at Mt. 
Abu has been observed by us over the last several years. Using a value ~ 15% for 
systematic error in the measured cosmic-rate is thus a fairly conservative estimate. 

The measured cosmic-ray rates have also been normalized to a zenith angle of 
0° and variation of the normalized rates as a function of run number is shown in 
Fig. 12c (NN-2 trigger ) and Fig. 12d (NN-3 trigger). In order to quantify departures 
from a flat shape, we have also calculated the by using the following expression 
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Fig. 12. Measured cosmic ray trigger rate as a function of zenith angle (a) data collected 
between November 10, 2005 - January 30, 2006 with NN-2 trigger configuration;(b) data 
collected between November 03, 2007 - March 8, 2008 with NN-3 trigger configuration. 
The spUne curve shown in theses figure is the Monte Carlo-based estimated cosmic-ray 
trigger rate. Normalized cosmic ray rate as a function of run number (c) November 2005 
- January 2006 data with NN-2 trigger (d) November 2007 - March 2008 data with NN-3 
trigger. The segmentation of the data into various spells are indicated by full vertical Unes. 
Each spell corresponds to one lunation period and the notation SP# stands for spell number. 



[33] 



o _ "f^" {RcR, - R^gr .... 



where nmn is the number of 10 minute runs, RcR,i is the normalized measured 
cosmic-ray rate with ARcR,i as the corresponding statistical error, fl'^^ is the rel- 
ative systematic uncertainty and R^^ is the normalized Monte Carlo simulated 
trigger rate. The resulting values of x^/dof ~ 605.4/597 (probability = 0.40) ob- 
tained for the NN-2 trigger observation spell and x^/do f ~ 662.6/611 (probability 
= 0.07) obtained for the NN-3 trigger observation spell clearly suggest that there 
is a reasonably good matching between the simulated trigger rates and the experi- 
mentally measured trigger rates. Presence of a systematic error at > 15% level for 
November 2007 - March 2008 data may also be contributing to the relatively higher 
value of for this data spell. 

The homogeneity of the triggers produced within the trigger field of view is another 
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test to check the performance of the trigger generation scheme. Although finding 
the frequency distribution of the triggered pixels would have been the best way to 
validate this, unavailability of the hit pattem data at the hardware level does not 
allow us to check this directly. To circumvent this, we determine the frequency dis- 
tribution of the image centroid coordinates by using the following approach. First, 
the data is corrected for inter-pixel gain variation and then subjected to the standard 
two-level image 'cleaning' procedure [27] with picture and boundary thresholds of 
6.5a and 3.0a respectively. In the clean Cherenkov image, if ith PMT with coordi- 
nates Xi, Hi registers a signal Sj, the coordinates of the image centroid are given by : 
< X >= J2i SixJ Y.i Si and < y >= J2i SiyJ Y.i Si- Considering the central ~ 2.25" 
X 2.25° trigger region and dividing it into 225 cells of size ~ 0.15° x 0.15° each, 
we bin the image centroid co-ordinates into these angular cells. Assuming that the 
image centroid is very close to the brightest pixel ( i.e. pixel with highest signal), 
the resulting centroid map can be faithfully used for checking the homogeneity of 
the triggers produced within the trigger field of view. Defining Nceii,i as number of 
events in the ith cell and < Nceii > as the average number of events in each cell 
we calculate percentage residual of events (ANceu,i) in each cell where ANceii.i 
={Nceii,i -< Nceii >)/< Nceii >. lu ordcr to ensure that images are robust and are 
contained within the camera we have also used the following cuts for selecting the 
events for this analysis :|< x >\< 1.2°, |< y >\< 1.2°, SIZE > ABOd.c where 
(6.5 digital counts =1.0 pe ) and number of pixels participating in the image > 4. 
Fig. 13 shows the results of the percentage residual events {ANceu,i) for the NN-2 
and NN-3 observation spells . A gray shade coding scheme has been used in this 
figure for representing the percentage residual events. The corresponding results 
obtained for the NN-3 trigger configuration by using data collected during Novem- 
ber 03, 2007 - March 8, 2008 is shown in Fig.l3b. In the NN-2 data ( Fig.l3a) we 
find that ~ 69% of the cells have | ANceia \ with in ~ 20%. The corresponding 
number of cells, for the NN-3 data ( Fig. 13b), is found to be ~ 65%. Presence of 
non-functional pixels in the camera (typically in the range of 3 to 7 in a total of 
225 pixels and whose location keeps on changing with time) as well as excessive 
trigger rate in some pixels due to presence of relatively bright stars are the 2 main 
causes because of which | ANceii,i \ is more than ~ 30% in some regions of the 
camera. 

Compatibility of the trigger rates with Poissonian statistics and homogeneity of the 
triggers produced within the trigger field of view clearly suggest that the overall 
performance of the telescope is quite stable and the Monte Carlo simulation esti- 
mates regarding effective collection area and 7-ray acceptance are also quite reli- 
able. The agreement between the expected and actual performance of the telescope 
has also been checked by comparing the expected and observed image parameter 
distributions and results of this study can be seen in [10]. The consistent detection 
of a steady signal from the Crab Nebula above ~1.2 TeV energy, at a sensitivity 
level of ~5.0 a in ~25 h, alongwith excellent matching of its energy spectrum with 
that obtained by other groups, also reassures that the performance of the TACTIC 
telescope is quite stable and reliable [27,29,31]. However, more work has to be 
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Fig. 13. (a) Percentage residual events ( ANceU,i ), for the NN-2 data, on the basis of the 
image centroid coordinates by following gray shade coding scheme, (b) Same as (a) except 
with NN-3 trigger configuration. The gray shade coding scheme followed for representing 
the percentage residual events is also given at the bottom of this figure. 

done with regard to relative calibration of the telescope ( along the lines suggested 
in [34]) so that unavoidable systematic errors as a result of temporal changes in the 
light collection efficiency, gain and, most importantly, atmospheric conditions can 
be accounted for in the analysis. 
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A detailed simulation study for the TACTIC telescope has been presented in this 
work for optimizing its trigger field of view and the trigger generation scheme. 
Different trigger regions with 49,81,121,169 and 240 pixels were considered in this 
work and the results suggest that a trigger field of 1 1 x 1 1 (~ 3.4" x 3.4°) is quite 
optimum for achieving maximum detection rate for 7-rays from a point source. 
With regard to optimization of topological trigger generation, it is found that using 
a NN-2 trigger methodology is a better option as against the NN-3 trigger config- 
uration when a source needs to be observed at zenith angles of > 25°. Reasonably 
good matching of the measured trigger rates (on the basis of ~ 207 hours of data 
collected with the telescope in NN-2 and NN-3 trigger configurations) with that 
obtained from simulations reassures that the procedure followed by us in estimat- 
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ing the threshold energy and detection rates is quite reliable. Furthermore, excellent 
matching of Crab Nebula differential spectrum with that obtained by other groups, 
reported by us in our earlier work, also reassures that the procedure followed by us 
for obtaining the energy spectrum of a 7-ray source is quite reliable. Attempts are 
also being made at present to improve the sensitivity of the telescope. Apart from 
using new CPCs ( modified CPC with square entry and a circular exit) for enhanc- 
ing the light collection efficiency of the camera, we are also trying to reduce the 
coincidence gate width. Both these modifications are expected to lower the thresh- 
old energy of the TACTIC telescope. To conclude, we firmly believe that there is 
considerable scope for a TACTIC like imaging telescope for monitoring AGNs on 
a long term basis. 
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